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Abstract. Let L be the Lie algebra of a simple algebraic group defined over 
a field F and let if be a split maximal toral subalgebra of L. Then L has a 
Chevalley basis with respect to H . If char(F) 7^ 2, 3, it is known how to find 
it. In this paper, we treat the remaining two characteristics. To this end, 
we present a few new methods, implemented in MAGMA, which vary from the 
computation of centralisers of one root space in another to the computation 
of a specific part of the Lie algebra of derivations of L. 



On the occasion of the distinguished birthdays 
of our distinguished colleagues John Cannon and Derek Holt 

1. Introduction 

1.1. The main result. For computational problems regarding a split reductive 
algebraic group G denned over a field F, it is often useful to calculate within its Lie 
algebra L over F. For instance, the conjugacy question for two split maximal tori in 
G can often be translated to a conjugacy question for two split Cartan subalgebras 
of L. Here, a Cartan subalgebra H of L is understood to be a maximal toral 
subalgebra, that is, it is commutative, left multiplication by each of its elements 
is semisimple (i.e., has a diagonal form with respect to a suitable basis over a 
large enough extension field of F), and it is maximal (with respect to inclusion) 
among subalgebras of L with these properties; it is called split (or ¥- split) if left 
multiplication by h, denoted ad^, has a diagonal form with respect to a suitable 
basis over F for every h £ H. Such a Cartan subalgebra is the Lie algebra of a split 
maximal torus in G. 

The conjugacy question mentioned above can be answered by finding Chevalley 
bases with respect to each split Cartan subalgebra, so the transformation from 
one basis to the other is an automorphism of L, and subsequently adjusting the 
automorphism with the normalizer of one Cartan subalgebra so as to obtain an 
element of G. In this light, it is of importance to have an algorithm finding a 
Chevalley basis (see Section LT731 for a precise definition). Such algorithms have been 
discussed for the case where the characteristic of the underlying field is distinct from 
2 and from 3. However, the latter two characteristics are the most important ones 
for finite simple groups arising from algebraic groups, so there is a need for dealing 
with these special cases as well. This is taken care of by the following theorem, 
which is the main result of this paper. 

Theorem 1. Let L be the Lie algebra of a split simple algebraic group of rank n 
defined over an effective field F. Suppose that H is an ¥ -split Cartan subalgebra 
of L. If L is given by means of a multiplication table with respect to an ¥-basis of 
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L and H is given by means of a spanning set, then there is a Las Vegas algorithm 
that finds a Chevalley basis of L with respect to H . IfF = ¥ q , this algorithm needs 
0~(n 10 log(g) 4 ) elementary operations. 

Here 0~(N) means 0(N(\og(N)) c ) for some constant c. Recall (e.g., from [P7] 
Introduction]) that arithmetic operations in F are understood to be addition, sub- 
traction, multiplication, division, and equality testing. If F is the field ¥ q of size 
q, these all take 0~(log(g)) elementary operations. Performing standard linear 
algebra arithmetic, that is, operations on matrices of size m, like multiplication, 
determinant, and kernel (solving linear equations), takes 0(m 3 ) arithmetic opera- 
tions. 

Better estimates than those of the theorem are conceivable, for instance because 
better bounds on matrix multiplication exist. However, our primary goal was to 
establish that the algorithm is polynomial in nlog(g). Moreover, in comparison to 
the dimension 0(n 2 ) of L or the estimate 0(n 6 ) for arithmetic operations needed 
for multiplying two elements of L, the high exponent of n in the timing looks more 
reasonable than may seem at first sight. 

The proof of Theorem [T] rests on Algorithm [TJ which is really an outline of an 
algorithm further specified in the course of the paper. The algorithm is implemented 
in Magma [3j . We intend to make the implementation public as a Magma package 
once the code has been cleaned up. 

The algorithm is mostly deterministic. However, in some instances where F is 
of characteristic 2 (such as Method [B 2 SC ] and the case where L is of type D4; see 
Sections 1531 13.51 and !3.6[) we use the Meat-axe (cf. [TTJ[T2]) for finding a particular 
submodule of a given module. For finite fields, the Meat-axe algorithm is analysed 
in [15] and [HJ Section 2] : irreducible submodulcs of a finite i-modulc of dimension 
m over F g can be found in Las Vegas time 0~(m 3 log(q)) (in Section ITT61 below it 
is explained why this result can be applied to Lie algebras). For infinite fields, 
Meat-axe procedures are known; however, we know of no proof of polynomiality in 
the literature. 

Algorithm [1] assumes that besides L and H the root datum R of the underlying 
group is known (see Section [L2|) . However, in Section [5] we show that this root 
datum can be determined by running the algorithm a small number of times. 

Thanks to the characterization of Lie algebras of split reductive algebraic groups 
described in Theorem [2] below, we can view the Lie algebras in Theorem [T] as 
Chevalley Lie algebras, which are defined below. 

Our treatment of Lie algebras and the corresponding algebraic groups rests on 
the theory developed mainly by Chevalley and available in the excellent books by 
Borel [T] , Humphreys [T3] , and Springer [TS] . Our computational set-up is as in [7] . 

1.2. Root Data. Split reductive algebraic groups are determined by their fields 
of definition and their root data. The latter is of importance to the corresponding 
Lie algebra and will therefore be discussed first. Throughout this paper we let 
R = (X, $, Y, $ v ) be a root datum of rank n as defined in [7]. This means X and Y 
are dual free Z- modules of dimension n with a bilinear pairing (■, •) : X x Y — > Z; 
furthermore, $ is a finite subset of X and <E> V a finite subset of Y, called the roots 
and coroots, respectively, and there is a one-to-one correspondence v : <& — > < f )V 
such that (a, a v ) = 2 for all a e 

If a £ $ then s a : x t— > x — (x, a v )a is a reflection on X leaving $ invariant and 
W = {s a I a S $) is a Coxeter group. Similarly, : y 1— > y — (a, y)a y is a reflection 
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on Y leaving <I> V invariant and the group generated by all these is isomorphic to 
W. In particular there are att,...,cei £ linearly independent in X <E> Q, such 

that <& = $+ U $~, where $+ = $ n (Nai H h Nai) and <f>~ = — <&+. The roots 

a±, . . . , oti and the coroots a\ , . . . , a ; v are called simple. The number I is called the 
semisimple rank of L (and of G). 

The pair (W, S 1 ), where S = {s ai , . . . , s ai }, is a Coxeter system. The Cartan 
matrix C of R is the ^ x I matrix whose entry is {ai,Oj). The matrix C is 
related to the Coxeter type of (W, S) as follows: s ai s aj has order where 

{ai,a)){otj,a{) 
4 

The Coxeter matrix is (mij)i<ij<i and the Coxeter diagram is a graph-theoretic 
representation thereof: it is a graph with vertex set {1, ...,£} whose edges are 
the pairs with > 2; such an edge is labeled rriij. The Cartan matrix 

C determines the Dynkin diagram (and vice versa). For, the Dynkin diagram 
is the Coxeter diagram with the following extra information about root lengths: 
(o!j,aJ) < (ctj,a.Y) if and only if the Coxeter diagram edge {i, j} (labelled my) is 
replaced by the directed edge (i, j) in the Dynkin diagram (so that the arrow head 
serves as a mnemonic for the inequality sign indicating that the root length of on 
is larger than the root length of ctj). 

A root datum is called irreducible if its Coxeter diagram is connected. A root da- 
tum is called semisimple if its rank is equal to its semisimple rank. Each semisimple 
root datum can be decomposed uniquely into irreducible root data. The Dynkin 
diagrams of irreducible root systems are well known, and described in Cartan's 
notation A„ (n > 1), B„ (n > 2), C n (n > 3), D„ (n > 4), E„ (n 6 {6,7,8}), F 4 , 
G2. The nodes are usually labeled as in [2]. 

For computations, we fix X = Y = Z™ and set (x, y) — xy T , which is an element 
of Z since x and y are row vectors. Now take A to be the integral / x n matrix 
containing the simple roots as row vectors; this matrix is called the root matrix of 
R. Similarly, let B be the I x n matrix containing the simple coroots; this matrix 
is called the coroot matrix of R. Then C = AB T and Z$ = ZA and Z$ v = TLB. 
For a G $ we define c a to be the Z-valued size / row vector satisfying a = c a A. 
In the greater part of this paper, including Theorems [T] and [2J we will let G be a 
simple group, so / = n. 




1.3. Chevalley Lie Algebras. Given a root datum R we consider the free Z- 
modulc 



L Z (R) = F©0ZA Q , 



where the X a are formal basis elements. The rank of L%{R) is n + |$|. We denote 
by [•,•] the alternating bilinear map L%(R) x Lz{R) — * L%(R) determined by the 
following rules: 

For y, z e Y : [y, z] = 0, (CBZ1) 

Forj/ey,a6$: [X a ,y] = (a,y}X a , (CBZ2) 

Forae$: [X_ a ,X a ] = a v , (CBZ3) 

For a, e *, a *±/* : ^] = { J^*** ^e' $ ' ( CBZ4 ) 



4 



AR.7EH M. COHEN AND DAN ROOZEMOND 



The N a ,p are integral structure constants chosen to be ±(p Qj( 3 + 1), where p a ,p is 
the biggest number such that a—p a ,pP is a root and the signs are chosen (once and 
for all) so as to satisfy the Jacobi identity. It is easily verified that N a p = — A_ Q ._£j 
and it is a well-known result (see for example [4j) that such a product exists. Lz{R) 
is called a Chevalley Lie algebra. 

A basis of Lz(R) that consists of a basis of Y and the formal elements X a and 
satisfies (CBZ1)-(CBZ4) is called a Chevalley basis of the Lie algebra Lz(-R) with 
respect to the split Cartan subalgebra Y and the root datum R. If no confusion is 
imminent we just call this a Chevalley basis of L%(R). 

For the remainder of this section, we let L%{R) be a Chevalley Lie algebra with 
root datum R, we fix X = Y = Z n , a basis of row vectors ei, . . . , e n of X, and a 
basis of row vectors /i, • ■ • ,/« of Y dual to ei, . . . , e n with respect to the pairing 
(•,•)■ Moreover, we let F be a field, we set hi = yt <g> 1, i = 1, . . . , rij and -ff = Y <g) F. 
Now tensoring Lz(R) with F yields a Lie algebra denoted Lf(R) over F, and the 
integral Chevalley basis relations (CBZ1)-(CBZ4) can be rephrased as: 



For i, j G {1, . . . ,n} : [hi,hj] 
For i G {l,...,n},a G $ : [X a ,/ii] 
For a £ $ : [A_ Q ,A Q ] 

For a,/j£f,a^ ±/3 : [* a ,-X/j] 

A Lie algebra is called split if it has a split Cartan subalgebra. The Cartan 
subalgebra H of each Chevalley Lie algebra L$(R) is split. The image of a Chevalley 
basis with respect to Y and R in Lf(R) is called a Chevalley basis of L with respect 
to H and R. The interest in Chevalley Lie algebras comes from the following result. 

Theorem 2 (Chevalley 5 ). Suppose that L is the Lie algebra of a split simple 
algebraic group G over F with root datum R — (A, $, Y, <I> V ), and that H is a 
split Cartan subalgebra of L. Then L = Lj?(R) and so it has a Chevalley basis 
with respect to H and R. Furthermore, any two split Cartan subalgebras of L are 
conjugate under G. Finally, if G is simple then R is irreducible. 

In light of this theorem, for the proof of Theorem [TJ it suffices to deal with 
Chevalley Lie algebras corresponding to an irreducible root datum. 

1.4. Some difficulties. So we will deal with the construction of a Chevalley basis 
for a Chevalley Lie algebra L over a field F, given only a split Cartan subalgebra 
H. Algorithms for finding such an H have been constructed by the first author 
and Murray [3] and, independently, Ryba [TB]. These algorithms work for char(F) 
distinct from 2 and 3, and partly for char(F) = 3. The first algorithm has been 
implemented in the Magma computer algebra system [3] . For now, we assume that 
we are also given the appropriate irreducible root datum R, but in Section [5] we 
argue that R can be found from L and H without much effort. The output of our 
algorithm is an ordered basis {X a , hi | a G 3>,i G {1, . . . , n}} of L$ (based on some 
ordering of the elements of $) satisfying (CB1)-(CB4). 

For fields of characteristic distinct from 2, 3, an algorithm for finding Chevalley 
bases given split Cartan subalgebras has been implemented in several computer 
algebra systems, for example Magma [3] and GAP [8]. For details, see for example 
[H Section 5.11]; the algorithm CanonicalGenerators described there produces 



0, (CB1) 

(a,fi)X a , (CB2) 

= E?=i(ei,a>i, (CB3) 

N a jjX a+l3 if a + e <J>, (QB4) 
otherwise. 



COMPUTING CHEVALLEY BASES IN SMALL CHARACTERISTICS 



5 



a Chevalley basis only up to scalars. The scaling, however, can be accomplished by 
straightforwardly solving linear equations. 

If, however, we consider Lie algebras of simple algebraic groups over a field F of 
characteristic 2 or 3, the current algorithms break down in several places. Firstly, 
the root spaces (joint eigenspaces) of the split Cartan subalgebra H acting on L are 
no longer necessarily one-dimensional. This means that we will have to take extra 
measures in order to identify which vectors in these root spaces are root elements. 
This problem will be dealt with in Section [3] Secondly, we can no longer always 
use root chains to compute Cartan integers (a,/3 v ), which are the most important 
piece of information for the root identification algorithm in the general case. We 
will deal with this problem in Section |4] Thirdly, when computing the Chevalley 
basis elements for non-simple roots, we cannot always obtain X a+ p from (CB4) 
by X a+ p = N [X a , Xp] as N a> p may be a multiple of char(F). This problem, 
however, is easily dealt with by using a different order of the roots, so we will not 
discuss this any further. 

1.5. Roots. Let p be zero or a prime and suppose for the remainder of this section 
that F is a field of characteristic p. We fix a root datum R — (X, &,Y, $ v ) and 
write L = Lf(R). We define roots and their multiplicities in L as follows. A root 
of H on L is the function 

n n n 

a : h ^ Ui)U, where h = 2J Ui ® U — khi, 

i—1 i—1 i—1 

for some a G here (a,yi) is interpreted in Z (if p = 0) or Z/pZ (if p ^ 0). Note 
that this implies that (a, h) := a(h) for h G H is completely determined by the 
values (a, yi}, i — 1, . . . , n. We write $(L, H) for the set of roots of H on L. 
For a G $(£, H) we define the root space corresponding to a to be 

n 

L-a = P| Ker(ad fti -a(hi)). 

i=l 

If a 7^ for all a G $ then L is a direct sum of H and its root spaces {L ff | a G $}. 
If on the other hand there exists an a G $ such that a = 0, then L is a direct sum 
of L = Cl(H) and a G $, a ^ 0}. 

Given a root a, we define the multiplicity of a in L to be the number of (3 G 
such that a = p. Observe that if a ^ the multiplicity of a G $(£, ff) is equal to 
dim(La). If a — this multiplicity is equal to dim(Lo) — n - Note that a 1— > a is a 
surjective map $ — > $(i, if), so in what follows we abbreviate $(£, if) to (f>. 

If each root has multiplicity 1, there is a bijection between $ and $. Our first 
order of business is to decide in which cases higher multiplicities occur. Observe 
that a = if and only if —a = so the multiplicity of the 0-root space is never 1. 
If char(F) = 2, then all nonzero multiplicities are at least 2 as a and —a coincide. 
Steinberg [HI Sections 5.1, 7.4] studied part of the classification of Chevalley Lie 
algebras L for which higher multiplicities occur (the simply connected case with 
Dynkin type A„, D„, Eg. 7,8 if char(F) = 2) in a search for all Lie algebras L 
with Aut(L/Z(L)) strictly larger than G. In Section [2] of this paper we prove the 
following proposition, which generalizes Steinberg's result to arbitrary root data. 
As the multiplicity of a root of H on the Lie algebra L of a central product of split 
reductive linear algebraic groups is equal to the minimum over all multiplicities of 
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ChevalleyBasis 

in: The Lie algebra L over a field F of a split reductive algebraic group, 

a split Cartan subalgebra H of L, and 

a root datum R = (X, $, Y, $ v ). 
out: A Chevalley basis B for L with respect to H and R. 
begin 

1 let = FindRootSpaces(£, H), 

2 let X = FindFrame(L, H, R, $, E), 

3 let i = IdentifyRoots(L, H, R, I>, X), 

4 let X°,H° = ScaleToBasis(L, H, R, X, i), 

5 return X°,H°. 
end 

Algorithm 1 . Finding a Chevalley Basis 

its restrictions to summands of the corresponding central sum decomposition of L, 
the study of multiplicities of roots can easily be reduced to the case where G is 
simple. 

Proposition 3. Let L be the Lie algebra of a split simple algebraic group over a 
field F of characteristic p with root datum R. Then the multiplicities of the roots 
in <E> are either all 1 or as indicated in Tabled 

In Table [T] the Dynkin type R of L and the characteristic p of F are indicated by 
R{p) in the first column. Further details regarding the table (such as the isogeny 
type of R appearing as a superscript on R{p)) are explained in the beginning of 
Section [2] This description uniquely determines the root datum R and hence the 
corresponding connected algebraic group G up to isomorphism; see [HI Chapter 9] . 

1.6. Computing ideals of Lie algebras. Finding an ideal / of a given Lie algebra 
L is equivalent to finding the submodulc / of the A- module L, where A is the 
associative subalgebra of End(L) generated by all ad x for x running over a basis of 
L. Hence, such an ideal / can be found by application of the Meat-axe algorithm 
to the A- module L. We will apply the Meat-axe only to modules of bounded 
dimension, so that the factor n 6 resulting from the occurrence of dim(L) 3 in the 
above-mentioned estimate for the Meat-axe running time when F = F ? plays no 
role in the asymptotic timing analysis. 

1.7. Algorithm Q3 In the remainder of this section we give a brief overview of the 
inner workings of Algorithm [T] It is assumed that L is isomorphic to L^(R). The 
FindRootSpaces algorithm consists of simultaneous diagonalization of L with 
respect to ad/^ , . . . , ad^ , where {hi, . . . , h n } is a basis of H. Its output is a basis 
E of iJ-eigenvectors of L and the set $ of roots of H on L. This is feasible over 
F because the elements are semisimple and H is split. As dim(L) = 0(n 2 ), these 
operations need time 0~(n 6 log(gr)) for each basis element of H, so the total cost 
is 0~(n 7 \og(q)) elementary operations. 

The algorithm called FindFrame is more involved, and solves the difficulties 
mentioned in Section [1. 41 bv various methods. The output X is a Chevalley frame, 
that is, a set of the form {FA Q | a G $}, where X a (a £ $) belong to a Chevalley 
basis of L with respect to H and R. If all multiplicities are 1 then FindFrame is 
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B 2 SC (2) 
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[Der] 
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[Der] 


F 4 (2) 
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[Der] 


G 2 (2) 


43 


[Der] 


B„ sc (2) ( n > 5) 


2™, 4(3) 


[C] 


all remaining(2) 


2 |»+| 


[A 2 ] 



Table 1. Multidimensional root spaces 



trivial, meaning that X = {Fx | x € E \ H} is the required result. The remaining 
cases are identified by Proposition^ and the algorithms for these cases are indicated 
by [A 2 ], [C], [Der], [B 2 SC ] in Tableland explained in Section El 

In IdentifyRoots we compute Cartan integers and use these to make the 
identification l between the root system $ of R and the Chevalley frame X computed 
previously. This identification is again made on a case-by-case basis depending on 
the root datum R. See Section [4] for details. 

The algorithm ends with ScaleToBasis where the vectors X a (a e $) be- 
longing to members of the Chevalley frame X are picked in such a way that 
X° = (X a ) a< z<s> is part of a Chevalley basis with respect to H and R, and a suitable 
basis H° = {hi, . . . , h n } of H is computed, so that they satisfy the Chevalley ba- 
sis multiplication rules. This step involves the solving of several systems of linear 
equations, similar to the procedure explained in [B], which takes time 0~(n 8 log(q)). 

Finally, in Section we finish the proof of Theorem [T] and discuss some further 
problems for which our algorithm may be of use. 



2. Multidimensional root spaces 

In this section we prove Proposition [3l but first we explain the notation in Table 
[T] As already mentioned, the first column contains the root datum R specified 
by means of the Dynkin type with a superscript for the isogeny type, as well as 
(between parentheses) the characteristic p. A root datum of type A3 can have 
any of three isogeny types: adjoint, simply connected, or an intermediate one, 
corresponding to the subgroup of order 1, 4, and 2 of its fundamental group Z/4Z, 
respectively. We denote the intermediate type by A3 . For computations we fix 
root and coroot matrices for each isomorphism class of root data, as indicated at 
the end of Section [L2l For A3, for example, the Cartan matrix is 
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As always, for the adjoint isogeny type A3 ad the root matrix A is equal to the 
identity matrix / and the coroot matrix B is equal to C . Similarly, for A3 SC we 

(2) 

have A = C and B = I. For the intermediate case A3 ; for instance, we take 



A = 












1 


I and B = 1 















It is straightforward to check that indeed det(A) = 2 = det(-B) and AB T = C. 

A root datum of type D„ has fundamental group isomorphic to Z/4Z if n is 
odd, and to (Z/2Z) 2 if n is even. The unique intermediate type in the odd case 
is denoted by D„ \ and the three possible intermediate types in the even case by 
d£\ Djr x \ and Dl n) . 

The multiplicities appear in the second column under Mults. Those shown in 
bold correspond to the root 0. For instance, for B2 SC (2) we have dim(CL(i/)) = 6, 
so the multiplicity equals 6 — 2 = 4. 

The third column, with header Soln, indicates the method chosen by our algo- 
rithm. Further details appear later, in Section [3J 

Assume the setting of Proposition [3l By Theorem [2] there is an irreducible root 
datum R = (X, $, Y, <1> V ) such that L = Lie(G) satisfies L = L W (R). Also, all split 
Cartan subalgebras H of L are conjugate under G, so the multiplicities of Lf(R) do 
not depend on the choice of H . For the proof of the proposition, there is no harm 
in identifying L with Lf(R) and H with the Lie algebra of a fixed split maximal 
torus of G. 

As all multiplicities are known to be 1 if char(F) = 0, we will assume that 
p := char(F) is a prime. We will write = for equality mod p (to prevent confusion 
we will sometimes add: mod p). We begin with two lemmas. 

Lemma 4. Let a, (3 £ <J>. Then a = (3 if and only if (c a — c l3 )A = 0. 

Proof. For h G H, by definition, (a, h) — (c a A,h) — c a Ah T . This implies that 
a = f3 if and only if c a Ah T = crAh for all h £ H, which is equivalent to 
(c Q - c^A = 0. □ 

Lemma 5. Let R\, R2 be irreducible root data of the same rank and with the same 
Cartan matrix C and denote their root matrices by A\ and A2, respectively. 

(i) //det(^42) strictly divides det(j4i), then the multiplicities in Lf[R±) are greater 

than or equal to those in ip(i?2)- 
(ii) If p /det(C), then the multiplicities of L^{R\) and Lp(i?2) are the same. 

Proof, (i). Without loss of generality, we identify the ambient lattices X and Y 
with Z™ and choose the same bilinear pairing (as in Section II. 2[) for each of the 
two root data R\ and i?2- The condition that det(^2) strictly divides det(Ai) then 
implies that the columns of A\ belong to the lattice spanned by the columns of A^. 
Hence A\ — A^M for a certain integral n x n matrix M. Thus (c a — cr )^2 = 
implies (c a — c^)A\ = (c a — c^AiM = 0, proving the lemma in view of Lemma[2] 

(ii). As det(C) ^ 0, the determinants of the coroot matrices B\ and B2 are nonzero 
modulo p, and A 1 = A 2 (B 2 B^ 1 ) and A 2 = A^B^ 1 ). It follows that (c a - 
c f 3 )A 2 = is equivalent to (c Q - c fJ )A 1 =0. □ 
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A typical case where part (i) of this lemma can be applied is when the adjoint 
and simply connected case have the same multiplicities, for then every intermediate 
type will have those multiplicities as well. It immediately follows from Lemma [S] 
that the root space dimensions are biggest in the simply connected case, and least 
in the adjoint case. Thus considering root data of the adjoint and simply connected 
isogeny types often suffices to understand the intermediate cases. Part (ii) indicates 
that in many cases even one isogeny type will do. 

The proof of Proposition [3] follows a division of cases according to the different 
Dynkin types of the root datum R. For each type, we need to determine when 
distinct roots a, (3 exist in $ such that a — [3. By Lemma[5]Jii), there are deviations 
from the adjoint case only if p divides det(C). 

As W embeds in Ng(H)/T, and acts equivariantly on $ and $ = $>(L,H), the 
multiplicity of a root a £ $ only depends on the VF-orbit of a £ $. By transitivity 
of the Weyl group on roots of the same length in it suffices to consider only 
a = ct\ in the cases where all roots in $ have the same length (A n , D„, Eg^s) and 
a = a.\ or a„ if there are multiple root lengths (B„, C„, F4, G2). 

In the adjoint cases, the simple roots ai, . . . , a n are the standard basis vectors 
ei, . . . , e„, since then the root matrix A and the coroot matrix B are / and C T , 
respectively. Similarly, in the simply connected cases, the simple roots a±, . . . ,a n 
are the rows of the Cartan matrix C, since then A = C and B = I. We write c = c 13 
so (3 = cA and either all c, e N or all Cj £ — N. 

We give the proofs of the cases where R is of type A„, B„, or G2. The other 
cases are proved in a similar way. For V a linear subspace of L and x £ L, we write 
Cy{x) for the null space of &d x on V, i.e., 

C v (ar) := {v £ V | [x,v] = 0}. 
2.1. A„(n > 1). The root datum of type A„ has Cartan matrix 



G 



/ 2 







-1 







2 -1 
-1 2/ 



and the roots are 

±(a j 



i}, k £ {j, ...,n}, 



■ + a k ), je{l,- 

where {ai, . . . , a n } are the simple roots, thus giving a total of 2 • hn(n + 1) roots. 

For the adjoint case, suppose aT = (3. Observe that all c, £ {0, ±1}. Since 
A = J, we must have ci = 1 and Cj = (J = 2, . . . , n), which implies either j) ^ 2, 
ci = 1, and C2 = • • ■ = Cn = 0, or p = 2, c\ = ±1, and C2 — ■ ■ ■ — c n = 0. Since we 
assumed (3 ^ ot\ we find p = 2 and /3 = —ol\, giving " £ n root spaces of dimension 
2. 

In the simply connected case the simple roots are equal to the rows of C, so that 
ai = j3 implies 2ci — C2 = 2, — c\ + 2c2 — C3 = —1, — Cj_2 + 2cj_i — Cj =0 for 
j = 4, . . . ,n, and — c„_i + 2c„ = 0. 

We distinguish three possibilities: c\ = 1, c\ = 0, and c\ = —1. If c\ = 1, then 
C2 = 0, so C2 = 0. As c±ai + - ■ ■+c n a n must be a root, this implies C3 = • • • = c n = 0, 
forcing (3 = of, a contradiction. 
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If ci = 0, then — C2 = 2, so that either p = 2 and C2 = 0, or p = 3 and c% = 1. 
In the first case, we find C3 = 1, giving a contradiction if n > 5 (because then 
C4 = and C5 = 1), a contradiction if n = 4 (because then the last relation becomes 
= — C3 + 2c4, which is not satisfied). Consequently, n = 3 and p — 2; the resulting 
case is discussed below. In the second case, where p = 3 and c 2 = 1, we find 
— I = 2 — C3, so that C3 = 0, giving a contradiction if n > 4 (because then C4 = I), a 
contradiction if n = 3 (because then the last relation becomes = — C2 + 2c3, which 
is not satisfied). It follows that n — 2 and p = 3; this case is also discussed below. 

If ci = —I, then —C2 = 4, so that either p = 2 and C2 = 0, or p = 3 and C2 = — f . 
In the first case, we find C3 = • • • = c n = 0, so /3 — —ot\. In the second case, we find 
that either n = 2 (the special case below), or C3 — 0, which leads to a contradiction 
if 7i > 4 (because then C3 = but C4 7^ 0), and also if n = 3 (because then the last 
equation becomes = — C2 + 2C3). 

We next determine the multiplicities in the two cases found to occur for A„ sc . 
For n = 3 and p — 2 we have 







A = C 



'o i cn 

1 1 I mod 2. 

,0 I 0; 



This gives al = 03, as well as a\ + ai — ct2 
ing for 3 root spaces of dimension 4. 
For n = 2 and p — 3 we have 



A = C = 



2 

-f 



-1 



■f «3 and «2 = cti 

I ^ I mod 3, 
Similarly, — cci = 



Q?2 



0:3, account- 



-a 2 



which implies a\ = 012 and aTj" = — (ai + 0^2). 
giving 2 root spaces of dimension 3. 

For the intermediate cases observe that by Lemma [SJi) we need only consider 
(n,p) — (2,3) and (3,2). But the former case has no intermediate isogeny types, 
and the latter case is readily checked to be as stated. This finishes the proof for 
the A„ case. 

2.2. B n (n > 2). The root datum of type B„ has Cartan matrix 



C 



/ 2 





Vo 



-1 







-2 
2/ 



and the roots are 

±(a 3 H h a{), 



±(atj H h + 

giving a total of 2 • \n(n 

In the adjoint case we have A 
■ = c, 



2ai + ■ ■ 

+ 1) + 2 



3 e {i, 

+ 2a n ), j G {1, 



,n}, 



,n},l e {j, ...,n}, 
,n-l},le O' + l, 

|n(n — 1) = 2n 2 roots. 

/. For the long roots, suppose oT = /?, so c\ = 1 
and C2 = ■ • • = c„ = 0. If ci = 1, then c 2 7^ (for otherwise /? = ai), which implies 
p = 2 and j3 = a\ + 2o/2 + • • ■ + 2a n . If ci = — 1, then p = 2, and either C2 = 0, 
which gives [3 = —a.\, or C2 7^ 0, which implies [3 — —a\ — 2a 2 — ■ ■ ■ — 2a n . In this 
case the long roots have multiplicities 4. 
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In the adjoint case, for the short roots, suppose = ft, so c n = 1 and c\ = ■ • • = 
c„_i = 0. This yields three possibilities for c„: If c n = —2, then p = 3, implying 
c n _i is either or —3, neither of which give rise to roots. If c n = — 1, then p = 2; 
now either c„_i = (yielding f3 = —a n ), or c n ~i = —2 (not giving any roots). If 
c n = 1 we must have c n _i = • • • = Ci = 0, giving the contradiction (3 = a n . This 
shows that p = 2 and all multiplicities are 2. 

In the simply connected case we have A = C. We will first consider n > 5, and 
then treat n = 2, 3, 4 separately. By Lemma EJii), we may assume p = 2. 

For the long roots, suppose ct[ = /3, so C2 = 0, ci + C3 = 1, and Cj_2 + Cj = 
(.7 = 4, . . . ,n). This forces C4 = 0. If ci = then c\ — and hence C2 = 0, so 
C3 = ±1. replacing /3 by j3 if needed, we may assume C3 = 1. As C4 = and C5 = 1, 
we must have C4 = 2 and c + 5 = 1, which is never satisfied by a root. If on the other 
hand ci = 1 then C3 = C4 = • • • = c n = 0, so f3 = —a% or [3 — ±(ai+2a2 + - • -+2a n ). 
This shows that, for n > 5, the multiplicities of j3 for /3 a long root are 4. 

For the short roots, suppose a^" = 0, so C2 = 0, Cj_2 + Cj =0 (J = 3, . . . , n — 1), 
and c n _2 + c„ = 1. If ci = 1 then C3 = 1, but since C2 = this contradicts that /3 is 
a root. If on the other hand c\ = 0, then C2 = C3 = ■ • • = c„_i = 0, so c„ = 1 and 
we find (3 — —a n . Hence, for n > 5, the multiplicities of /3 for j3 a short root are 2. 



If n — 2 then 

C 



2 -2\ = / 
-1 2 " 1 



If ~a\ = (3 we have ci = 0. Since —2 < C2 < 2 we must have either C2 = (hence 
(3 = — ai), or C2 = ±2 (hence ci = ±1), giving (3 = ±ai or (3 = ±(ai + 20:2). If 
on the other hand 02 = /? we find C2 = 1 hence (3 — ±ct2 or (3 = ±(qi + a^)- This 
shows that £?2 SC has 2 root spaces of dimension 4 if p = 2. 



If n = 3 then 



/ 2 -1 ' 
C= -1 2 -2 
\0 -1 2, 





1 




) ■ : 





• 




1 





From a straightforward case distinction on the roots of B3 and the fact that ai = 
03 we immediately see that aT — S3 = a\ + 2a2 + 2«3, 02 = ol\ + oli + a.3 = 
(X2 + 2oi3, and ai + 012 = «2 + «3 = «i + «2 + 2«3. This gives the 3 required root 
spaces of dimension 6. 



If n = 4 then 

/0 1 0\ 
10 10 
10 0' 
\0 1 0/ 

From a straightforward case distinction on the roots of B4 and the fact that 
a! = 03, we find aTf = S3 = 0:3 + 2a4 = ai + 2^2 + 2a3 + 2a4, as well as 02 = 
Q!i + «2 + us = a.\ + o?2 + «3 + 2a4 = CK2 + 2«3 + 2«4 and ai + a2 = Q!2 + 013 = 
a2 + 2a3 + 2«4 = ai + «2 + 2a3 + 2a4. The remaining 32 — 24 = 8 roots (±(oy + 
• • ' + a n), j = 1, . . . , 4) are in 2-dimensional spaces, giving 2 4 , 8 3 , as required. 



C = 



( 2 


-10 \ 


-1 


2-10 





-1 2 -2 


\o 


0-12/ 
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2.3. G2- The root datum of type G2 has Cartan matrix 




and the roots are 

±ai, ±(ai + as), ±(2ai + a-z), (6 short roots) 
±02, ±(3ai + 0:2), ±(3ai 4- 2a2), (6 long roots) 

giving a total of 12 roots. As det(C) = 1, we take A — I. All components of c 
are in {—3, . . . , 3}, so all components of the differences ot\ — (3 and 012 — (3 are in 
{—4, . . . , 4}. Hence, if multidimensional root spaces occur, we must have p < 3. 

If p = 3 we see 3ai + — 02 = — (3oti + 2^2) and — (3ai + 02) = — «2 = 
3ai + 2^2, and the remaining 6 roots all have distinct root spaces. 

If p = 2 we find ol\ + a2 = 3cci + 012, oT = 3ai + 2«2 and 02 = 2«i + a%, giving 
3 root spaces of dimension 4. 

This finishes the proof of Proposition [3] 

3. Finding Frames 

Let L be a Chevalley Lie algebra over F with root datum R, a fixed split Cartan 
subalgebra H, and given decomposition E into root spaces with respect to the 
set <& = $(I/,iJ) of roots of H on L. In this section we discuss the procedure of 
Algorithm [T] referred to as FindFrame. It determines the set X = [FX a | a G $}, 
i.e., the one-dimensional root spaces with respect to $, to which we refer as the 
Chevalley frame. Note that we do not yet identify the root spaces: finding a suitable 
bijection between $ and the Chevalley frame X is discussed in the next section. 
We set p = char(F). 

We require that R be given, since we execute different algorithms depending on 
R, for example B2 ad needs [C] whereas B2 SC needs [B2 SC ]. 

For p = 2, we use the procedure described in Section [3. II to find the frame once 
we have computed all spaces ¥X a + ¥X^ a for a £ $. To this algorithm we will 
refer as [A2]. As an auxiliary result, this procedure stores the unordered pairs 
{{a,— a} I a G to be used in the IdentifyRoots procedure discussed in 

Section 2] (notably, the proof of Lemma [5]) . 

The general method in characteristic 2 is to reduce the root spaces of dimension 
greater than 2 to such 2-dimensional spaces, and apply [A2]. For this purpose, and 
for the two cases of characteristic 3, we distinguish three general methods: 

• [C]: Given two root spaces M,M' compute Cm(M') to break down M. 
Often, but not always, dim(M') = 2. An example of this method is given 
in Section [3~2l 

• [Der]: Compute the Lie algebra Der(L) of derivations of L, and calculate in 
there. This is a useful approach if Der(L) is strictly larger than L, for then 
we can often extend H to a larger split Cartan subalgebra, so we find new 
semisimple elements acting on the root spaces. Examples of this method 
are given in Sections 13.31 and 13.41 

• [B2 SC ]: The case where R(p) = B2 SC (2) is slightly more involved than the 
other cases because a — for some a G <&. We use the Meat-axe to split 
the action of the long roots on the short roots. Examples of this method 
are given in Sections 13.51 and 13.61 
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The method chosen depends on the root datum R and the characteristic p, as 
indicated in the third column of Table [TJ 

3.1. A2 in characteristic 2. First, we consider the Lie algebras L with R(p) = 
A2(2), as this procedure is used inside various other cases. The isogeny type of the 
root datum is of no importance here. For clarity, we write a, (3 for the two simple 
roots of the root system of type A 2 . 

As indicated in Table [TJ we have 3 root spaces of dimension 2. They correspond 
to (Xy, A_ 7 )f for 7 G {a, (3, a + j3}. Without loss of generality we consider La — 
{X a , X- a }f and L-g — (Xp, X^p)^. Observe that the squared adjoint action ad^ Q 
of X a sends any element of L-s to zero; [A a , [-X^-Xg]] — [A a , JV a ^X a -\-p] — as 
2a + (3 £ $, and [X a ,X-p] = since a - (3 <&. Similarly, ad^_ Q {Irs) = 0. 

However, the quadratic action ad^ of a general element x = t\X a + £ 2 X_ Q 
(ti,t2 G F, both nonzero) of L„ does not centralise L-g. Indeed: 

[x,[x,X p ]] = htiQX-a^XcXpW + lX^lX-^Xp]]) 
= tit 2 N- ata+l 3N a ^Xp, 

which is nonzero since N— a a +/3 and N a ,p are both equal to 1 modulo 2. 

Recall that we are given L„ and L-s. Fix a basis r\ , r 2 of La and consider the 
element x = r% + tr 2 , where t G F. It follows from the above observations that 
ad^(L^g) = if and only if a; is a scalar multiple of X a or X^ a , so in order to find 
the frame elements among the Fa; for t G F we have to solve 

= [x,[x,y]] = [ri + tr 2 ,[ri+tr2,y]] 

= [r x , [ri,y]]+t ([r x , [r 2 ,y]\ + [r 2 , [r u y]}) + t 2 [r 2 , [r 2 ,y]], 

for every y G L-a in the unknown t. We know there is a solution as H is split. Solving 
this system is equivalent to solving a system of 2 • 3 = 6 quadratic equations in t 
(note that the [ri, [rj, y]} are in (L-q)l, which is at most 3-dimensional) . If F = ¥ q , 
solving such a quadratic equation is equivalent to solving log(q) equations in log(g) 
variables over F 2 (as p = 2 is fixed), requiring 0~(log(g) 3 ) arithmetic operations, 
or 0~ (log(g) 4 ) elementary operations. 

For more general Lie algebras L, the solutions for Lie subalgebras of type A 2 
normalized by H will be part of a Chevalley frame. These parts can be found 
inside any two-dimensional root space V G E provided there is at least one other 
two-dimensional root space V' G E such that (V, V') l is of type A 2 . So, if all root 
spaces in E are 2-dimensional and F = ¥ q , this method needs 0(n 2 ) root spaces 
V to be analysed (at a cost of 0~(n s log(g) 4 ) each), so that X will be found in 
0~(n 10 log(g) 4 ) elementary operations. 

3.2. G2 in characteristic 3. Secondly, we consider the Lie algebra L — Lf(G 2 ) 
of the root datum of type G 2 over a field F of characteristic 3. By Proposition [3] 
there are 8 root spaces. It is readily verified that dim(L ff ) = 1 if a is a short root 
and dim(La) — 3 if a is a long root of <&. In particular, the short root spaces belong 
to X and it remains to split the two long root spaces. 

Consider one of the two three-dimensional root spaces in E, say V = VX a2 + 
¥X 3ai +a 2 +^X^ 3ai -2a 2 ■ The left multiplications on V by the short roots are easily 
obtained from (CB1)-(CB4); these are given in Table [21 

Although we have not yet identified the roots, we can identify the three pairs 
of one-dimensional root spaces {FX Q ,FX_ Q }, for a G $ short, since L— is the 
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Xa 2 


^3ai+a2 


X. 


3ai 


-2a 2 


X ai 


Xai +ct2 


























X ai -\-ci2 








X- 


-2ai 


— Q2 


V 

— ai — a2 


—X ai 












X2 ai +a 2 












ai 


X—2ai -Q2 





—X ai 










Table 2. Part of the G2 multiplication table 



unique one-dimensional root space with root —a. From this observation and Table 
[2] it follows that we can obtain the triple FXp (f3 £ {«2, 3ai + ct2, — 3ai + 202}) as 
follows: 

¥X a2 = C v (L 2ai+a2 +L_ 2ai _ a2 ) 7 
VX 3ai+a2 = G v (L ai+a2 +L_ ai _ a2 ), 

For the other three-dimensional space, the same approach is used. This completes 
the search for the Chevalley frame X. 

3.3. D4 in characteristic 2. Thirdly, we consider the Lie algebras with Dynkin 
diagram of type D4 over a field F of characteristic 2. As mentioned in Section [2l 
there are three cases: 

L ad : the adjoint root datum (12 two-dimensional root spaces), 

L sc : the simply connected root datum (3 eight-dimensional root spaces), 

Lv 1 ), L( 3 \ l/ 4 ); the intermediate root data (6 four-dimensional root spaces). 

The three intermediate root data all give rise to the same Lie algebra up to iso- 
morphism (by triality), so we will restrict ourselves to the study of L ad , L sc , and 
LW. It is straightforward to verify that L ad has a 26-dimensional ideal I ad , linearly 
spanned by X a (a € $), {a( +ot% + a\) ® 1, and a 2 <£> 1. This ideal can be found, 
for example, by use of the Meat-axe. 

Similarly, L sc has a 2-dimensional ideal / (spanned by (a^ + a( ) <S> 1 and (a^ + 
at) ® !)• Let I sc — L sc /I be the 26-dimensional Lie algebra obtained by computing 
in L sc modulo /. Finally, L^ 1 ' has a 1-dimensional ideal / (spanned by 014 ® 1), and 
a 27-dimensional ideal /' (spanned by 014 (£> 1 and X a , a S $). We let 1^ = I' /I. 
Again, the 26-dimensional ideal is easily found by means of the Meat-axe. 

Thus we have constructed three 26-dimensional Lie algebras: I ad , I sc , and 
By results of Chevalley (cf. [TH Part 2, Cor. 2.7]) they are isomorphic, so from 
now on we let I be one of these 26-dimensional Lie algebras. The Lie algebra / is 
simple. Its derivation algebra Der(I) is a Lie algebra of type F4, and thus has 12 
two-dimensional root spaces and 3 eight-dimensional root spaces. 

Using a procedure similar to the one for G2 over characteristic 3 described in 
Section ^. 21 we can break up the eight-dimensional spaces of E into two-dimensional 
spaces, giving us 24 two-dimensional spaces. These two-dimensional spaces may 
then be broken up into one-dimensional spaces by the procedure [A2] . The last step 
in the process is "pulling back" the relevant one-dimensional spaces from Der(J) to 
I. But this is straightforward, since I is an ideal of Der(7) by construction. 
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3.4. G2 in characteristic 2. As noted in [19j Section 2.6], in the exceptional case 
R(p) = G2(2), the Lie algebra L is isomorphic to the unique 14-dimensional ideal 
of the Chevalley Lie algebra L A of adjoint type A3 over F. 

In particular, Der(L) contains a copy of L A . Wc use this fact by finding a 
split Cartan subalgebra H' inside Crj er (£)(£f) so that H C H' . For then we can 
calculate the Chevalley frame X A inside the Lie subalgebra {L, if')Der(i) of Der(L) 
with respect to H' , which is of type A3 by the above observation. 

The Chevalley frame X of L is now simply the part of X A that lies inside L. 

3.5. B2 SC in characteristic 2. We consider the Chevalley Lie algebra L of type 
B2 SC over a field F of characteristic 2 with split Cartan subalgebra H = ¥hi + ¥h 2 - 
This is a particularly difficult the automorphism group of L is quite big: 
Aut(L) = G x (F+) 4 [lQl Theorem 14.1], where G is the Chevalley group of adjoint 
type B 2 over F and F + refers to the additive group of F. As a consequence, there 
is more choice in finding the frame than in the previous cases. 

To begin, we take L to be the (0, 0)-root space of H on L, and L\ to be the (1, 0)- 
root space of H on L. It is easily verified that Lq — (H, X± ai , X±f ai+2 a 2 ))N (that 
is, the linear span of H and the long root elements) and L\ = (X± a2 , A±( Qi+Q2 ))f 
(the linear span of the short root elements). We proceed in three steps. 

[B2 sc .l]. The subalgebra Lq has Dynkin type Ai © Ai. We may split it (non- 
uniqucly) into two subalgcbras of type Ai using a direct sum decomposition pro- 
cedure. This is a procedure that can be carried out with standard linear algebra 
arithmetic for a fixed dimension (6, in this case); see e.g., [9] Section 1.15]. 

[B2 SC .2]. Let A be one of these subalgebras of Lq of type Ai. Assume for the sake 
of reasoning that A = (X± ai )L, the Lie subalgebra of L generated by X ai and 
X^ ai . Since [A, L{\ = L\ we may view L± as a four-dimensional A-module, and 
hence apply the Meat-axe [TTJ Q~2] to find a proper irreducible A-submodule M of 
L\. This will be a submodule of the form 

M = (t\X a . 2 + t 2 X 

We take 61,62 t° be a basis of M, and add Ga^) and C^(6i) to X. These 
two spaces are indeed one-dimensional and coincide with the original ¥X± ai if 
b\ e ¥{tiX a2 + i2A_ Ql _ a2 ) and 62 G ¥(tiX ai+a2 + tiX- a2 ). This exhibits part of 
the freedom of choice induced by the factor (F + ) 4 in Aut(L). 

We repeat this procedure for both subalgebras of type Ai found in the first step. 
The result is the part of the Chevalley frame X inside Lq- In fact, due to our 
method, we can make an identification of the long roots ±ai, ±(ai + 2012) with 
the four elements of X found. In what follows we will work with such a choice so 
that we have the elements ¥X ai , FA_ ai , FX ttl+ 2a 2 , FA_ Ql _2 Q2 in X as well as 
the correspondence with the roots in $ suggested by the subscripts. 

[B2 SC .3]. We find the part of X inside L\ as follows. FX ai+a2 coincides with 
CL 1 (FA Ql ,FA ai +2a 2 )- Having computed this element of X 7 we finish by taking 

FA Q2 = [FA Q , 1+Ct2 , FX—aiU 

¥X- ai - a2 = [FA Q2 , FA_ Ql _2a 2 ], 

FA_ Q2 = [FA Ql _ Q2 ,FA ai ]. 

This completes the search for X in the case B2 SC (2) and establishes that its running 
time is 0~(log(g)). 
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3.6. C n sc in characteristic 2. We consider the Chevalley Lie algebra L of type 
C„ sc over a field F of characteristic 2. Here n > 3, so that the multiplicity of is 
strictly larger than 4. Let h z be a basis of the 1-dimensional centre of L, inside 
the split Cartan subalgebra H of L. This case is a generalisation of the B2 SC case 
described in Section [3~5l We again take Lq to be the 0-root space of H on L, so that 
Lq is 3n-dimensional and consists of H and the root spaces corresponding to the long 
roots. Similar to the previous case, Lq = Ai © • • • © Aj (n constituents), and again 
the decomposition is not unique. We describe how to find such a decomposition. 

We let T be the set of Q) four-dimensional root spaces (cf. Table [1]). In the root 
system of type C„ each of these corresponds to the four roots ±ej ± ej for some 
i,j G {1, . . • , n} with i ^ j. Our first task is to split Lq into subalgebras of type 
Ai in a way compatible with T . To this end, we let T be the graph with vertex set 
T, and edges / ~ g whenever / ^ g and [/, g] ^ 0. 

Let A be a maximal coclique of L of size n — 1, so that A consists of n — 1 
elements of J- such that [/, g] = for all /, g € A. This means that, for a particular 
i G {1, . . . , n}, the set A C T corresponds to those four-spaces in T that arise from 
the roots ±ej ± 6j, where j € {1, . . . , n} \ {i}. Let A = L — A, so that A contains 
precisely the four-dimensional spaces corresponding to ±ek ± Q with fc, I ^ i. 

Now compute the centralizer A in Lq of all spaces in A. Then A coincides with 
(X± 7 ,7 V ® l,h z )f for the long root 7 = 2ei. Using a direct sum decomposition 
procedure we find the Lie subalgebra A' of A such that A — A' © ¥h z , where 
A' = (X± 7 ,7 V <g> The subalgebra A' is one of the type Ai constituents of Lq 
we are after. Thus, by repeating this procedure for each maximal coclique of T of 
size n — 1, we obtain a decomposition of Lq into n subalgebras of type Ai. We will 
denote by A the set of these n subalgebras. 

Now we continue as in the B2 SC case: For each element of A we use the procedure 
labelled [B2 SC .2] to find suitable elements FX±~ f for X. For each four-dimensional 
space K £ T we then use distinct S\,Si € A satisfying [A, Si] ^ 0, [A, 5 2 ] 7^ and 
these FA± 7 to execute a [B2 SC .3] procedure. Thus, we find the part of the frame 
inside A. 

If n = 3 splitting Lq has to be done in a slightly different way, but as this is only a 
slight modification of the algorithm we will not go into details here. This completes 
the Chevalley frame finding in the case C„ sc (2). Its running time involves 0(n 2 ) 
executions of parts of the algorithm of Section [3. 51 which is however dominated by 
the time 0~(n 10 log(g) 4 ) needed for method [A2]. 

We summarise the results of this section. 

Proposition 6. Given L, H, R, the set $ 0/ roots of H on L, and the root spaces 
E, the Las Vegas procedure FindFrame finds a Chevalley frame. For ¥ = ¥ q , it 
runs in time 0~(n 10 log(g) 4 ). 

Proof. As mentioned in Section 11.71 this procedure is trivial in all cases except 
those mentioned in Table [TJ and for each of the cases in Table Q] we have presented 
a solution. Recall that |$| < dim(L) = 0(n 2 ). 

The timing of Method [A 2 ] is dealt with in Section lBTTl which produces the bound 
stated in the proposition. 

Method [C] concerns 0(n 2 ) instances of standard linear algebra arithmetic on 
spaces of bounded dimension, and so its running time is dominated again by time 
spent on the [A 2 ] method. 
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Method [Der] involves the computation of parts of the Lie algebra of derivations. 
Computing the full Lie algebra of derivations in instances like D n sc (2) would take 
running time 0~ (n 12 log(g)). However, we only carry out this procedure for Lie 
algebras of bounded dimension (the bound being 28, which occurs for type D4) 
or compute the part of Der(L) that leaves invariant H and the corresponding de- 
composition into root spaces (which reduces the running time to 0~(n 8 log(g))). 
Therefore, the stated bound suffices. 

Finally, according to Table [TJ Method [B2 SC ] with unbounded n only occurs in 
the cases treated in Section [3. 61 where the time analysis is already given. □ 



In this section we clarify Step 3 of the ChevalleyBasis algorithm [T] The 
routine IdentifyRoots takes as input a Chevalley Lie algebra L, a split Cartan 
subalgebra H of L, the root datum R and the set of roots 3> = §(L, H), and the 
Chevalley frame X found in the previous step (Section [3]). It returns a bijection 
1 : $ — > X so that, up to scaling, (X a ) a6 $ will be the root element part of a 
Chevalley basis. 

An important tool to make this identification are the Cartan integers (a,/3 v ). 
Cartan integers may be computed using root chains; see, for instance, [4]. 

Lemma 7. Let a,/J £ $. Suppose p and q are the largest non-negative integers 
such that a — p(3 6 $ and a + q[3 £E 4>. Then (a, /? v ) = p — q. 

We may use this lemma by computing such a chain in the set of roots <f> cor- 
responding to the Chevalley frame X = {¥X a | a S $}. However, as these roots 
are computed from the Lie algebra L over F itself, they live in the n-dimensional 
vector space F" rather than over Z n . 

A straightforward verification of cases for Chevalley Lie algebras arising from 
root systems of rank 2 shows that the chain can simply be computed in terms of 
the roots over F™, except if the characteristic is 2 or 3. So in the latter two cases, 
we a different method for computing (a,/3 v ) is needed. 

Lemma 8. Suppose that L = Lf(R) is a Chevalley Lie algebra with respect to a root 
datum R — (X, $, Y, $ v ) over the field F of characteristic 2 or 3. Let H be the stan- 
dard split Cartan subalgebra of L. Suppose furthermore that X a , X- a , Xp, X-p are 
four vectors spanning root spaces corresponding to a, —a, (3, —(3 6 $, respectively, 
and a 7^ ±/?. 

// <J> is simply laced, then (a,/3 v ) = P — Q, where 



( if[X_ ,X a ] = O 

P = { 1 tf[X-f,,X a ]?0,[X-f,,[X-(,,X a ]]=0 

I 2 if[X^,[X_ ,X a ]]^O 

f if[Xp,X a } = 

Q = < 1 if[Xp,X a ]^O i [Xp,[X p ,X a ]]=0 

{ 2 if[X p ,[Xp,X a ]]^0 



4. Root Identification 




If $ is doubly laced and char(F) ^ 2, then (a,{3 v ) 



P — Q, where 
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Proof. For any 7,7' G let p 7 , 7 ' and 97,7' be the biggest non-negative integers 
such that 7 — P7,7'7' G $ and 7 + q ln >^' G Recall from (CB4) and [4] that, if 
7 + 7' G $, then [Xy,Xy/] = A 7i7 /A 7+ y, where A 7i y = ±(p ln > + 1). 

If $ is simply laced, the subsystem of $ generated by ±a, ±/3 is of type A1A1 or 
of type A 2 . Then a + /3 G $ implies a — /3 G" so A Qi/3 = ±1. Similarly, Np <a = ±1. 
This means that, regardless of the characteristic, we can reconstruct p a ,/3 and g Qj( 3 
by the procedure described in the lemma, and thus compute (a,/3 v ) = p a ,p ~ Qa,/3 
by Lemma [7] 

If $ is doubly laced and char(F) 7^ 2, the subsystem of <J> generated by ±a, ±/3 is 
of type A1A1, A2, or B2. (Note that G2 never occurs inside a bigger Lie algebra.) 
In the first two cases the previous argument applies, so assume ±a, ±/3 generate 
a subsystem of $ of type B2. Similarly to the previous case, if a + (3 G $ then 
a — 2/3 G" $, so that N a p,Np iQ G {±1,±2}. In particular, since char(F) 7^ 2, we 
find that both N a p and Np >a are nonzero, so that we can reconstruct p a ,p and g Qj( 3 
by the procedure described in the theorem, and thus compute (a, /3 V ) = p Q „8 — la, p 
by Lemma [7] □ 

Lemma 9. Suppose that L is a Chevalley Lie algebra over F with respect to a root 
datum R = (A", $, Y, $ v ), H is a split Cartan subalgebra of L, and X a and Xp are 
two root elements whose roots with respect to H are a and (3 for certain a,/3 G 
Suppose, furthermore, that one of the following statements holds. 

(1) char(F)£{2,3}, 

(2) $ is simply laced, 

(3) $ is doubly laced and char(F) 7^ 2. 

Then (a,/3 v ) can be computed from the available data in 0~(n 10 log(g)) elementary 
operations. 

Proof. Observe first of all that the case where a = (3 is easily caught, for example 
by computing dim((FA Q , FJ^f). Obviously then (a,/3 v ) = 2. 

Moreover, we can distinguish the case where a = —(3 as follows. If char(F) 7^ 2 
we may simply test whether a = —(3. If on the other hand char(F) = 2, we find 
the sets {{7, —7} | 7 G ( I >+ } as an auxiliary result of the algorithm FindFrames 
described in introduction of Section ETTl If a = —(3, then of course (a,/3 v ) = —2. 

So assume a 7^ ±(3. Now if (UJ holds we compute (a, (3 V ) from the roots a and 
(3 using Lemma UJ as mentioned earlier. 

Suppose, therefore, @ or ([3]) holds. We can find FA_ Q and FA_^ either simply 
by considering {7 | 7 G $} (if char(F) 7^ 2) or as an auxiliary result of FindFrames 
(if char(F) = 2). This leaves us in a position where we may apply Lemma [8l and 
thus find (a,/3 v ). 

Finally, the time needed does not exceed the time needed for standard linear 
algebra arithmetic for each pair of roots, that is, (9~(n 4 • n 6 \og(q)). □ 

The last lemma enables us to compute Cartan integers in many cases. For 
the cases not covered by Lemma |H] we proceed as follows to construct a direct 
identification 1. 

• B„(2): The short root spaces generate an ideal, I say, of L found by the 
Meat-axe, and the root eigenspaces of H that do not lie in / belong to long 
roots. These root spaces generate a subalgebra of type D„. This Lie algebra 
is simply laced, so the root identification problem can be solved there. This 
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identifies the long root spaces. Now, for i = 1, . . . , n, let the short root 7, 
be cti + ct.;+i + ■ ■ - + a n and let «o = or + 2a2 + 2«3 + • • • + 2a n be the (long) 
highest root. Observe then that [A Qo , A_ 7l ] = A 72 and [X ao ,X^ l2 ] — A 7l , 
and A_ 7l and A_ 72 are the only short root elements that do not commute 
with X ao . This fact, together with the set of pairs {{7,-7} | 7 <E < I )+ } 
obtained in FindFrames, allows us to find X± 11 and X± l2 . Note that we 
have to execute this procedure at most twice, since there are only elements 
of X that could be identified with A_ 7l , and the other short root elements 
are fixed once X_ 7l is fixed. 

The other short root elements may now simply be found by using rela- 
tions such as [X 7i , A_ Q J = X li+1 . 

• C„(2): The short root spaces generate an ideal of L of type D„, so we 
execute a similar procedure as in the previous case. 

• F4(2): The short roots generate generate an ideal of L of dimension 26 
which together with the Cartan subalgebra H gives a 28-dimensional sub- 
algebra of type D4, allowing the same procedure as before. 

• 62(3): Similarly to the previous cases, we use the fact that the short roots 
generate an ideal of L of type A2, which is again simply laced. 

• 62(2): As described in Section [3^1 the manner in which the root spaces in 
L A correspond to those in L is completely determined. Therefore, we may 
use the roots identified in L A , which is simply laced, to identify the roots 
in L. 

These methods lead to the following conclusion. 

Proposition 10. Given L over F, H, R = (X, $,Y, $ v ), the set $ of roots of 
H on L, and a Chevalley frame X , the routine IdentifyRoots finds a bijection 
1 : 4> — > X such that for all a,/?£$, a ^ ±/3, 



For F = F q , the routine needs 0~ (n 10 log(q)) elementary operations. 

Proof. Lemma |H] shows that in many cases we can compute Cartan integers. To 
this end, we need to compute (a,/3 v ) for all 0(n 4 ) pairs of roots, and every com- 
putation of this type involves at most 6 multiplications in L, requiring a total 
of 0~(n 4+6 log(q)) elementary operations. Once these numbers are computed, it 
takes 0(n 4 ) steps to select a set of simple roots and subsequently to complete the 
bijection between $ and X. These last two steps use techniques similar to those 
described by De Graaf [5J Section 5.11]: the creation of a set of simple roots n 
starts with taking an arbitrary root to be the first member of IT. We then itera- 
tively pick a suitable additional simple root f3 having Cartan integer (a,/3 v ) < 
with the members a of IT. This proves that we can make the required bijection in 
0~(n 10 log(g)) time for the cases covered by Lemma [HI 

For the remainder of the proof, we can restrict ourselves to the cases not covered 
by Lemma [S] Here the procedure described provides 1 directly, so we only need 
prove the last assertion of the proposition. As G2(2) is directly reduced to a case 
already treated, it needs no further consideration. In each of the remaining cases, 
we need to compute a subalgebra or an ideal of L. Although this is hard in general, 
the fact that we have already found the Chevalley frame X and the fact that the 
subalgebra or ideal is a sum of elements from X imply that the computations take 




t(a + (3) if a + (3 e $ and N a ,p ^ (mod p), 
{0} otherwise. 
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Table 3. Algorithm [T] timings 



0~(n w log(g)) elementary operations. A bijection i! from the relevant subsystem 
of <& to the subset of X of root spaces lying in the ideal may then be identified in 
time 0~ (n 10 \og(q)) . Finally, extending i! to the entirety of $ is a straightforward 
task, requiring only standard linear algebra arithmetic in L. 

This shows that we can make the required bijection in time stated for all cases. 

□ 

5. Conclusion 

As discussed in Section 11.71 the more difficult steps of Algorithm [1] are FlND- 
Frame and IdentifyRoots. In Sections [3] (Proposition!!]) andl4l fProposition flOl) 
we established that these steps can be dealt with in time 0~ (n 10 log(<7) 4 ). This 
proves Theorem [T] for a given root datum. 

We emphasize that this estimate is only asymptotic. Additionally, in Table [3] we 
present timings of Algorithm [T] for various root data and for four different fields: Q, 
GF(17), GF(3 3 ), and GF(2 6 ). The times given are in seconds, for the most time- 
consuming root datum with the specified Lie type. Input for the algorithm were 
a Chevalley Lie algebra and its splitting Cartan subalgebra, to which a random 
basis transformation was applied, and the root datum. The timings were produced 
using Magma 2.15-5 on an Intel Core 2 Quad CPU running at 2.4 GHz with 8GB 
of memory available, although only one core and at most 2.7GB of memory were 
used. 

As hinted at earlier, Algorithm [T] can easily be used to produce an algorithm 
that takes only L and H and produces the root datum R and a Chevalley basis. 
To see this, note first that, because H is given and the underlying algebraic group 
is assumed to be simple, we may use dim(_ff) = rk(i?), the dimension of L, and the 
classification of simple Lie algebras to narrow down the root system to one or two 
possibilities (or three, but only if dim(L) = 78 and dim(_ff) = 6). Therefore, the 
number of possible root systems it at most 3. 
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Second, given a root system, the number of possible root data is small as well. 
If the root system is not of type A or D, the number of possible isogeny types is at 
most 2. If the root system is of type D the number of possible isogeny types is at 
most 5, as explained in Section[2] So suppose <& is of type A„, and fix p — char(F). 
Note that the fundamental group is Z/(n + 1)Z. Since two root data for A„ lead to 
isomorphic Lie algebras if both have the same exponent of p in [X : 1A>\ , we need 
consider at most log p (n+ 1) + 1 = 0(log(n)) different isogeny types. Thus, in order 
to identify the correct root datum, we run Algorithm [T] a sufficiently small number 
of times for the polynomiality bound given in the theorem to remain intact. This 
finishes the proof of Theorem [TJ 

A problem hinted at, but not solved satisfactorily, is finding a split Cartan 
subalgebra of a Chevalley Lie algebra L if p = 2. Nevertheless, verifying that a 
given subalgebra is indeed a split Cartan subalgebra is easy. So our results are 
still useful, since one is often able to obtain such a subalgebra by other means, for 
example as part of the original problem. Moreover, experimental implementations 
of randomized algorithms for finding split Cartan subalgebras look promising. We 
intend to publish about these algorithms in forthcoming work. 

A primary goal in writing the Chevalley basis algorithm is to use it for conjugacy 
questions in simple algebraic groups G or finite groups G(¥ q ) of rational points over 
W q . One of the complications in this application is the fact that the group Aut(£) 
may be much larger than G(¥ q ). For this purpose, a method is needed to write 
an arbitrary automorphism of I as a product of an element from G(¥ q ) and a 
particular coset representative of G(¥ q ) in Aut(L). Such a method is in [7] and is 
also used in [6] . 

Once Algorithm [T] completes succesfully we have a certificate for a Lie algebra 
to be of type R: when presented with a candidate Chevalley basis X°, H°, we only 
need to carry out the straightforward and quick task of verifying that X°,H Q is 
indeed a Chevalley basis for L with respect to H and R. In this way, our work also 
contributes to a recognition procedure for modular simple Lie algebras. Obviously, 
Algorithm Q] can be used for establishing an isomorphism between two Chevalley 
Lie algebras over the same field and of the same root datum. 
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